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Abstract. For the adhesive contact of elastomers, we propose expressions to quantify 
the impact of viscoelastic response on effective adhesion energy as a function of contact 
edge velocity. The expressions we propose are simple analytical functionals of the creep 
response and should be suitable for experimental data analysis in terms of measured 
rheologies. We also emphasize the role of the coupling between local stress field 
at the contact edge and the macroscopic remote loading (far field). We show that 
the contrast between growing and receding contacts originates from the impact of 
viscoelastic response on coupling, while the separation process at the contact edge is 
similarly affected by viscoclasticity in both cases. 



1. Introduction 

Technological applications of soft polymers are multifarious. An especially noteworthy 
property is their high compliance resulting in remarkable contact behaviours. This is 
why they have been used for about a century for tyre manufacturing, for example. 
As a result the adhesive contact of elastomers has been thoroughly explored over the 
last decades. The problem neatly combines mechanical response of the materials and 
interfacial properties P, El [3l HI [5], [6]. Over the years a very thorough theoretical 
understanding has been reached, to such an extent that adhesive elastomeric contacts are 
routinely used for the assessment of surface properties: the contact will be established 
and the force necessary to pull out the sphere monitored as a function of the (decreasing) 
contact radius. By applying the JKR model, an evaluation of the adhesion energy [71[8], 
characteristic of the interface is obtained. However, beyond the JKR model, the adhesive 
contact of elastomers hides some more complex features. 

1.1. Receding contacts and enhancement of adhesion 

While he JKR model for elastic bodies involves the reversible adhesion energy (denoted 
w here), it is usually found that for the rupture of elastomeric adhesive contacts, the 
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measured adhesion energies are considerably larger than w. This enhancement results 
from some viscoelastic dissipation which occurs in the polymeric material but originates 
from the tensile cohesive stresses involved in the reversible adhesive process [HI [TO] . 
The question is the quantitative evaluation of this dissipative contribution given the 
viscoelastic response of the elastomer. 

In this respect, the elastomeric adhesive contact problem, and the closely connected 
viscoelastic crack problem, have received a lot of attention over the years [5l [11], [121 [131 
[m [iSl [iSl [13, [TBI [ini [20]. To describe the viscoelastic dissipation at the edge of the 
contact, the standard technique is to introduce a cohesive zone (Fig. [1]) to include the 
reversible adhesive stresses as boundary conditions to the viscoelasticity problem. The 
argument is well-known: for a linear viscoelastic system, an infinitely sharp tip would 
lead to infinite strain rate - and no dissipation - whatever the crack tip velocity |21j . 
In contrast, introducing the range of the cohesive interactions, there appears a zone 
with finite size (a-c in Fig [Tj) over which the interaction stresses do work. Then the 
viscoelastic material experiences a finite strain rate resulting in velocity dependent 
dissipation through the viscoelastic response. 

The loose statement that the viscoelastic crack response scales with the dissipative 
component of the linear mechanical response is often found [6l[22]. A more quantitative 
result can be derived from the formulation by Schapery and Greenwood [HI [151 
However this formulation has yet only been applied to idealized materials. Indeed in 
practice the numerical calculations are complex [231 CZ] 5 so that approximate schemes 
have been developed [121 [131 HH [21] • Another theory which takes into account the 
dissipative processes generated by a cohesive zone has been proposed [251 [26]. Heralding 
the fact that there is no such direct relation with the dissipative function, it has become 
known under the denomination of the "viscoelastic trumpet". In the present paper we 
propose a new approximate expression accounting for the adhesion energy enhancement 
by viscoelastic dissipation in a receding contact and demonstrate that the form is simple 
to use and more accurate than some standard approximations. 

1.2. Growing contacts and suppression of adhesion 

An even more striking manifestation of viscoelastic dissipation in the adhesive contact 
of elastomers is the contrast between growing and receding contacts. As just mentioned 
above, for contact rupture (a receding contact), an increased adhesion energy is recorded. 
However, during loading, for an adhesive contact growing fast, the contact morphology is 
similar to an adhesionless (Hertz) contact (Fig.[2]c) and the equilibrium JKR-like contact 
(Fig. [2] a) is very slowly reached when the penetration, or the load, is held constant: for 
an elastomer which equilibrates instantaneously in a bulk mechanical test, it typically 
requires ten minutes for an adhesive contact to reach equilibrium [27]. This latter effect 
is less known because the suppression of adhesion is somehow less appreciated than its 
enhancement. 

However, the loss of adhesion signals another interesting question. Indeed the 
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cohesive zone process determines the local stress field at the the edge of the contact 
(or crack tip); the latter in turn must be matched with the far field, a process we call 
coupling (Fig. [1]). As we will show in this paper, the issue of the coupling in the adhesive 
contact of elastomers is simple for receding contacts, but not for growing contacts. For a 
growing contact, it was first quahtatively discussed by Johnson and Greenwood [Hj, then 
more explicitly dealt with by Schapery [15] at various degrees of approximation. These 
results where subsequently adopted by Johnson and Greenwood [28l[T7]. In the present 
paper, we explicitly exhibit simple approximate expressions for the coupling between 
local contact edge stress field and remote, far field, loading. For a growing contact, the 
velocity dependence of the coupling demonstrates how the viscoelastic process quells 
the adhesion by smothering the coupling. 

In brief, the emerging picture is that the measured adhesion energy is an effective 
adhesion energy (or apparent adhesion energy or apparent work of adhesion) and that 
this effective adhesion energy (denoted Ge// here) may be suppressed (growing contact) 
or enhanced (receding contacts) as a result of viscoelastic dissipation. 

1.3. Strategy - Effective adhesion energy, viscoelastic materials and elastomers 

For an elastomer, a soft polymeric material which does not flow (rubber-like material), 
the mechanical response is controlled by the molecular chain entropy [29]. At longer 
experimental time scales, elastomers exhibit a non-dissipative response. At smaller 
time scales, the viscoelastic behavior associated with the polymeric glass transition 
appears. However, as the glass transition temperature of elastomeric materials is low 
(typically below -30 °C), and because of the time-temperature superposition principle, 
the viscoelastic time scale is far much shorter than the second at room temperature. 
As a result the rationale behind the use of the JKR model for elastomers is as follows: 
the bulk of the system effectively behaves as an elastic solid and assumes the response 
of the relaxed material, while viscoelastic effects act at the local scale, near the contact 
edge, inducing dissipation [28] . 

This picture gives a a flrm ground to the notion of effective adhesion energy. Let 
us consider two geometrically identical systems, one ideally elastic, and the other one 
elastomeric; and, in addition, that the long time modulus of the elastomeric material is 
equal to the modulus of the perfectly elastic material. Due to dissipation of the elastomer 
at the contact edge, identical far flelds (same penetration, same applied force) require 
different reversible adhesion energies for the elastic system: a larger for the receding 
contact, a smaller for the growing contact. Conversely the effective adhesion of the 
elastomeric system is equal to the reversible adhesion energy which is necessary to 
produce identical far fields in the elastic system. 

In the present paper we revisit the theory of the adhesive contact of elastomers from 
a different viewpoint. Indeed there exists another set of relevant contact theories, which 
deal with the fully viscoelastic systems, in which the contact timescale is of the order 
of the stress relaxation time, so that the process of stress relaxation in the bulk must 
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Figure 1. Schematics of tiie macroscopic contact and the cohesive zone along with 
the couphng between the far field on the one hand, and the energy transfer within the 
cohesive zone on the other hand. 



be explicited. By restricting such theories to elastomeric behaviour - the elastic bulk 
approximation- a consistent picture of viscoelastic adhesive contacts may be elaborated. 

For the full viscoelastic contact problem, the solution proposed by Falsafi [30] was 
later criticized by Hui who put forth a more consistent theory [21] for both growing 
and receding contacts. Other solutions were later proposed based on the powerful 
Hankel transform technique [3ll [321 [331 [SI [351 [3B] pioneered by Sneddon [37] . Indeed, 
some years ago, using the Sneddon method, we have proposed another theory for the 
viscoelastic adhesive contact [38l [39l [40]. Two main results were obtained: 1) an 
analytical expression for the viscoelastic crack response; 2) a general relation between 
local crack tip stress field and far field relevant for the coupling problem in a fully 
viscoelastic system. 

It is this theory which is the starting point for the present investigation of 
elastomeric adhesive contacts where we propose a new formulation for the effective 
adhesive response of elastomeric materials. For both growing and receding contacts, 
expressions are identified for the viscoelastic crack response (Eq. [TT!) and for the effective 
adhesion (Eqs. [T9] and [2T1) . These expressions are cast in the form of simple analytical 
functionals of the creep function (Eqs. [H [9] and [H]). We expect they will be useful 
in the analysis of experimentally determined effective adhesion energies based on the 
knowledge of the creep function. The relations which these expressions bear to other 
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Figure 2. Schematics of tlie adhesive contact for an elastomer.5jKR is the elastic 
JKR penetration for an adhesion energy w. For small velocities, the system is close 
to equilibrium and receding and growing contacts are close to the elastic (JKR) limit 
with the relaxed modulus. For large velocities, strong deviations connected to the 
viscoelastic dissipative processes appear, which are quantified in the effective adhesion 
energy. 



elastomeric contact theories are also discussed. 
2. Model 

2.1. Elastic Adhesive Contact 

We first briefly review the use of cohesive zones in adhesive contact problems with purely 
elastic systems. 

2.2. Cohesive zone - Elastic material 

In the cohesive zone the work of the tensile stresses during crack propagation induces a 
transfer between adhesion and mechanical energy [JT]. It has been shown that there is 
little impact of the details of the cohesive zone on the crack tip behaviour [42] so that the 
choice of a given cohesive zone model is not expected to signiflcantly affect the results. 
The cohesive zone model may be chosen for convenience and while Greenwood [17] used 
Maugis' [13] we prefer Greenwood's own "double Hertz" stress distribution [H] (Fig. [1]). 
Integrated over the cohesive zone size eo = c — a, the work of the cohesive stresses 
(~ (To) in the surface displacements (~ ao/E*) can be calculated self-consistently. 
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Here the elastic parameters are the Young's modulus E and the Poisson ratio u, while 
E* = £'/(! — z/^) is the reduced elastic modulus. It is equal to the adhesion energy 
w [ISl m] , so that for the double- Hertz model [S] 

^ = ^^1? (elastic) (1) 
2.3. Local stress field distribution 

The stress field at the edge of the contact zone (but inside, r < a, Fig. [1]) is controlled by 
the cohesive zone stress field (a < r < c). The resulting enhancement of the local tensile 
stress field is proportional to the stress intensity factor K. In the Sneddon formalism 
the form which appears naturally for the stress intensity factor is [38] : 

g{a) - --(ToV^ae^ (2) 



while K = —2g{a) I ^[Txa. 

The energy transfer between the physical process of adhesion and the local contact 
edge stress field (Eq. [1]) can be written 

w = (elastic) (3) 

which we identify with the familiar expression for the energy release rate 

w = ^ (4) 
2E* ^ ^ 

2.4- Coupling with the far field 

In the elastic case, the coupling between the far field and the contact edge stress field 
is given by 

2 

SjKR- ^Hia) = -f;-9{a) (elastic) (5) 

where 5h{o) is the adhesionless contact radius. The coupling coefficient between contact 
edge stress field and far field is C = 2/ E*. The force is 

FjKR{a) = Fnia) + Aag{a) = E^Fh^o) + 2aE\5jKR - 5/f(a)) (6) 

where Fh = E*Fh is the adhesionless contact force. In anticipation of section [3l Fh is 
introduced because it contains the geometrical dependence of the (adhesionless) contact 
force and leaves out the mechanical response E*. Such a separation is possible for 
homogeneous half-spaces only [3l]. Inserting Eq. [3] the familiar JKR equations are 
recovered [1]. 
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2.5. Far field - flat punch displacement 

More precisely Eqs. [5] and [6] allow for an interpretation in terms of an additional fiat 
punch displacement Sfp{a) = 5jkr{o) — 6H{a) which produces the stress field singularity 
proportional to g{a) at the crack tip [H Hi]. Let us start from an adhesionless contact 
{g{a) = 0): the penetration 5 is a function ^/^(a) of the contact radius a. Then we turn 
adhesion on and assume some adhesive process subsumed in g{a) < 0. Eq. [5] expresses 
the fact that to keep the contact radius constant while turning the adhesion on, one 
needs to pull back the sphere by a fiat punch displacement Sjkb^ — Snia) proportional 
to g{a). The resulting additional force term in Eq. [6] is the fiat punch force (the fiat 
punch displacement times the fiat punch stiffness 2aE*). 

Alternatively, the same interpretation can be formulated as follows: when adhesion 
is turned on for a given penetration 5 jkr, the contact radius increases so that Eq. [5] 
is fulfilled. The cohesive stresses tend to extend the contact zone in a manner similar 
to wetting. This tendency for the cohesive stresses to expand the contact zone clearly 
announces a possible role of creep in the coupling between local crack tip stress field 
and far field in the case of viscoelastic materials (section 13.21) . 

3. Viscoelastic Adhesive Contact 

For a linear viscoelastic material [15], following Schapery [TT] and Greenwood [17] in the 
usual restrictions on the definition of the response for contact models, we introduce the 
creep function (j){t) [38] which provides the delayed deformation in response to the history 
of applied stresses. Denoting Eoo the relaxed modulus (for the long time response) and 
Eq the instantaneous modulus (the short time response), we have 0(0) = 2/Eq* and 
0(+oo) = 2 /Eoo*. k is the contrast between relaxed and instantaneous moduli. An 
elastomeric material has k = E^/ Eq << 1. 

3.1. Cohesive zone model, creep and effective crack tip compliance 

For a fast crack, the bonding or debonding processes will be effectively elastic and the 
modulus involved will be the instantaneous modulus. Similar considerations apply to 
the slow crack, where the relevant modulus will be the relaxed modulus. The latter is 
smaller by a factor k resulting, for identical cohesive stresses, in larger deformations and 
proportionally smaller cohesive zone size for a slow crack (Eq. [3]) More generally, for 
a viscoelastic system, the cohesive zone size e depends upon crack tip velocity V and 
for intermediate velocities, the crack behavior will be characterized by some effective 

§ Our intuition tends to correlate sharp necking on the crack edge with smaU cohesive zones. This is 
an error: the necking is a resuh of the far field while the cohesive zone size is relevant locally at the 
crack tip. Indeed the cohesive zone size e is directly proportional to the square of the stress intensity 
factor g{a) (Eq.[TO|). Note also that for an elastomeric contact, g{a) only indirectly drives the flat punch 
displacement 5 — Snia), as will be shown below (Eqs . [T3l and fT5|) . As a result for similar cohesive zone 
sizes, a strong (Fig. [2]b) or nearly absent (Fig.[2]c) necking may be recorded at the macroscopic scale. 
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crack tip compliance of the material (denoted 0i here) and an intermediate value of the 
cohesive zone size. This is what a viscoelastic crack model must quantify. 

Let us define the dwell time tr (or transit time) as the time it takes for the 
moving crack to cross its own imprint e = c — a (Fig. [1]). In the quasi static crack 
tip approximation we assume the tip shape is convected, i.e. the crack velocity V does 
not change appreciably during a time comparable with the dwell time tr- We simply 
have V = etj. and we want to calculate e{V), or rather e{tr). 

In the viscoelastic case, we have shown [38j that the crack tip processes can be 
described by the following effective crack tip compliance: 

w = —o-o^e{tr)(pi{tr) (viscoelastic) (7) 
8 

where (piifr) is either (pi^opifr) or (pi^difr) with 

<Pi,op{t) = Jq^^ ^ r)(f){T)dT (opening crack) (8) 



and 



0i,d(i) = 77 / T(l){T)dT (closing crack) (9) 
t'^ Jo 

Locally the contact edge stress field, as characterized by the stress intensity factor 
g{a) is [38] 



TT 



g{a{t),tr) - --ao^2a{t)e{tr) (10) 
As a result we have 

w = 01 (^r) (viscoelastic) (11) 

na 

where 0i(tr) = 4>i,op{t) or (pi^dit) for an opening or closing crack. 

For practical purposes Eqs. [H]and[n]for the effective crack tip compliance 0i(tr-) may 
well be the main result of the paper. However, to be of use, it must first be connected to 
the measurable parameter, the effective adhesion, a task we complete only in section 1X51 

Here we will discuss some representative quantitative results for the effective crack 
tip compliances Eqs. [H] and O We can visualize the effective crack tip compliance by its 
impact on the cohesive zone size, which is directly (Eqs. [1] and [7j) 

ho 01 (tr) 

We consider the special case where the creep function is a three-element model 



( Appendix A ), which is frequently used as a coarse rendering of viscoelastic behavior. 
T is the creep time and the ratio of long time to instantaneous elastic modulus is 
k = 0.01. The normalized cohesive zone sizes for opening or closing cracks are plotted 
on Fig. [3l The results are expressed as a function of the "fast crack" characteristic 
velocity V = eo/T. As expected Eqs. E] and [9] predict a transition between the relaxed 
cohesive zone at low velocity and the instantaneous cohesive zone at high velocity. The 
transition occurs around a normalized velocity of 1. The minor offset between the 
opening and the closing cracks is simply due to the different cohesive stress history 
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Figure 3. Cohesive zone size as a function of velocity for opening and closing cracks. 
Creep follows an exponential time dependence (three-element model) with k = 0.01. 
The viscoelastic cohesive zone processes are nearly identical for closing or opening 
cracks. 
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Figure 4. Dwell time as a function of velocity for opening and closing cracks. Creep 
follows an exponential time dependence (three-element model) with k = 0.01. 



experienced by the cracks moving in opposite directions. Altogether however the two 
cases are very similar and the curves collapse when a scaling by a factor of about 2 is 
applied on the velocity. 

3.2. Coupling between local and far field 

3.2.1. Penetration For a growing contact, we found that (Eq. 14 in |38j ) 

5{t) — duiait)) ~ (l)o(tr)g{a{t),tr) (growing contact) (13) 

where 

Mt) = l fdT<Pit-T) (14) 
t Jo 

Of course, since the contact is growing, it is the expression for a closing crack (Eq.lH]) 
which must be taken for (pi in Eq. [TTJ 
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At this level of approximation, for a growing crack, the coupling coefficient between 
crack tip stress field and far field is C = 0o(^r)- The coupling coefficient directly refiects 
the creep acting over the cohesive zone during time tr and its impact on the penetration. 
As the effective crack tip compliance, the coupling coefficient exhibits a strong variation 
(Fig. [5]) with dwell time tr (or conversely with crack velocity. Fig. Hj). U tr is small - 
or the contact edge velocity large -, creep is inefficient and C ^ 2/E*q: the additional 
flat punch displacement is smaller by a factor k than the flat punch displacement for 
a relaxed material (Fig. [2] c) the necking characteristic of the JKR adhesive contact 
is nearly absent and contact growth appears to proceed as a virtually non-adhesive 
contact. If tr is large - or contact edge velocity small -, C ~ 2/E*oo and the additional 
flat punch displacement is large: we revert to the completely elastic, relaxed, adhesive 
contact (Fig. [2] a). 

For a receding contact 

2 

S{t) — Sniait)) — — — g{a{t),tr) (receding contact) (15) 

where 0i is for an opening crack (Eq. [8]) since the contact is receding. Eq. [15] turns out 
to be the equation for the elastic case (Eq. [5]) with the relaxed modulus. The coupling 
coefficient between crack tip stress field and far field is C = 2/E*oo, which is independent 
from tr (Fig. [5]). 

For a receding crack, creep plays no role in the coupling. The surface displacements 
resulting from creep occur outside of the contact zone and move away from it so that 
they will not affect the far field. The coupling between far field and contact edge 
stress field is as in the usual JKR model (Fig. [5]). Of course creep still plays a role in 
the (independent) self-consistent crack tip problem (section \2.2\i . which determines the 
effective crack tip compliance. 

3.2.2. Force To complete the picture we also study the force. Expressions for the force 
are more complex because they result from the spatial integral of the normal surface 
stress. In both cases, since the bulk is fully relaxed, the first term on the RHS is the 
Hertz force with the relaxed modulus. The additional flat punch term is the product of 
the (relaxed ) contact stiffness with the actual flat punch displacements 5{t) — Sniait)) 
as given by Eqs. [13] and [15] respectively. 
For a growing contact: 

Fit) ^ E^^Fniait)) + ^ct>o{tr)Aag{a{t)) (16) 

Note that the flat punch force term is controlled by the speciflc coupling coefficient for 
the growing contact. 

For a receding contact: 

F{t) ~ E^'^FHiait) + Aa{t)g{a{t)) (17) 

which is the force equation for the elastic case, with the relaxed modulus. 

II It is the flat punch displacement for an elastic system with the instantaneous modulus. Of course, 
the force (Eq. I16p is not the force for the instantaneous modulus. 
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3.3. Effective adhesion 

We now consider the effective adhesion Geff as defined in section 11.21 It can be 
calculated comparing Eqs. [I3] and [I5] with Eq. [5], or Eqs. [I6] and [I7| with El In both 
cases g{a) is determined by Eqs. [11] or [3] as appropriate. 

For the receding contact (using the opening crack expression for the effective crack 
tip compliance Eq. [8]) the effective adhesion is simply the normalized cohesive size 



so that 



w er 



^eff _ (Pl,op 



oo 



:i8) 



(receding contact) (19) 



W 4>l,op{tr) 

where 0i^op(oo) = 2/E*oo- For contact growth (with the closing crack expression for 
the effective crack tip compliance Eq. [H]), due to the non trivial coupling coefficient 
C = (poitr) 

-^^^^^%i{tr)- (20) 



w \ 2 / e 



Coo 



so that 



(growing contact) (21) 



W (f)l,cl{oo)(j)i^cl{tr) 

where 0i,d(oo) = 

Typical results are plotted on Fig. [6] for a three-element model with k = 0.01. 
Despite nearly identical cohesive zone responses for a given velocity (Fig. [3]), the effective 
adhesion energies for growing and receding contacts separate at higher velocities to reach 
their respective limit values equal to k and 1/k respectively. This contrast is due to 
the strong decrease of the coupling coefficient for growing contacts at higher velocities 

(Fig. ED. 



4. Discussion 



With Eqs. [19] and [21] we have provided expressions for the effective adhesion of 
elastomeric contacts for receding and for growing contacts. These expressions are very 
general and can easily be applied to arbitrary creep functions. We now compare our 
results to predictions and approximations of the literature. 



4.1. Receding contact and effective adhesion energy 

Based on earlier work by Schapery [12], Greenwood and Johnson [T3] have proposed 
a simple approximation. In this creep function approximation the effective crack tip 
compliance is equal to the creep function evaluated at some effective dwell time t*. By 
assumption the effective dwell time t* is proportional to the dwell time t,., that is t* = ptr 
for an opening crack. The approximation is completely defined if the proportionality 
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Figure 5. Normalized coupling coefficient C between crack tip stress field and far field 
as a function of contact zone velocity for growing and for receding contacts. Creep 
follows an exponential time dependence (three-element model) with = = 0.01. 




Figure 6. Effective adhesion energy for growing and for receding contacts. The nearly 
identical crack tip processes (Fig. [S]) result in very different effective adhesion energies 
because of the dependence of the far field/local field coupling coefficient with contact 
edge velocity (Fig. [5]). The result of the product rule is also shown (see [Appendix B[ ). 



coefficient p is specified. A value of 1/3 lias been proposed [mill]. In 2004 Greenwood 
performed numerical calculations to improve the evaluation of the effective dwell time 
t*. He proposed p ~ 1/4 for a three-element viscoelastic solid with k = 0.02 and a 
Dugdale-Maugis cohesive zone [43] . 

In the same paper [T7] he also showed that the creep function approximation is 
only valid at shorter dwell times (or high contact edge velocity): at longer dwell times 
(or low contact edge velocity), the creep function approximation for the effective crack 
tip compliance is inaccurate and an effective dwell time can no longer be defined. 

Here we have proposed a different type of expression for the effective crack tip 
compliances (Eqs. [8] and [9]). Compared to Greenwood's creep function approximation, 
which is based on the creep function itself, but evaluated at an effective dwell time, our 
expressions emerge from higher order developments and assume the form of moments 
of the creep function. For the three-element model, our effective crack tip compliances 
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Figure 7. Normalized effective compliances (pi for opening and closing cracks as 
predicted by Eqs.[5]and[ni), and comparison with the creep function approximation |14| 
(fc = 0.01). Our results match the creep function approximation at shorter dwell times 
(high contact edge velocities), where it is known to be valid, and also captures the 
deviations evidenced at longer dwell times (slower velocitites) known from the earlier 
numerical calculations [17| . The coupling constant (po as predicted by Eq. [14] and its 
short time expansion are also displayed. 



can be explicited analytically (Sec. Appendix A ). To compare with the creep function 



approximation, 0i op is plotted on Fig. [3 We find that this function tracks the creep 
function at short dwell times (or high contact edge velocities), allowing for explicit 
evaluation of the effective dwell time. Expansion of Eq. IA.4I shows that p = 1/3. This 
result is (accidentally) identical to Schapery's [12], and similar to Greenwood's |17j . 
Indeed explicit values may differ due to different cohesive zone types [I9] and different 
models for the viscoelastic response. 

However what also clearly appears (Fig. [7j) is that for longer dwell times (or lower 
contact edge velocities) the effective crack tip compliances converge more slowly to 
their asymptote than the creep function. This result is in full agreement with the trend 
identified by Greenwood using extensive numerical calculations [17] and demonstrates 
the improved validity of the higher order expansion. We believe that our approximate 
expressions retain most of the simplicity of the creep function approximation but provide 
better evaluations of the effective response especially at low contact edge velocities. 

4.2. Growing contact and effective adhesion energy 

It is apparent that the viscoelastic processes at the contact edge are very similar for 
receding and for growing cracks (Fig. [31 Indeed, in terms of effective adhesion energy, our 
results (Eqs. [T3]and [T5l) demonstrate that the difference between receding and growing 
contacts (Fig. [6]) lies in the coupling (Fig. [5]), not directly the contact edge viscoelastic 



crack response (Appendix B). At higher velocities, nearly identical viscoelastic crack 
responses for growing or receding contacts result in very contrasted effective adhesions 
because the coupling coefficients between local stress field and far field (Fig. [5]) differ 
by a factor k depending on the direction of contact edge motion. Similar results may 
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be found in the literature, especially in [13] (Eqs. 63) and [T7] (par. 3 and appendix) 
but we believe that stronger emphasis on the coupling, such as conveyed by Eq. [21] and 
Fig. [5], best brings out the workings of the growing contact. 



5. Conclusion 



We have proposed approximate analytical expressions (Eqs. [19] and [2T]) for the effective 
adhesion energy of elastomers. These expressions are easy to use because they involve 
simple functionals (Eqs. [HI M and [Hj) of the creep function. 

We have shown that our expressions match the existing approximate model, the 
creep function approximation, for high contact edge velocities, and improve on it at low 
contact edge velocities. In this regime, our model correctly predicts the deviation of the 
compliance from the creep function approximation which has been observed numerically. 

We also propose a consistent form for the effective adhesion energy of growing 
contacts, where emphasis is laid on the role of the coupling between local stress field 
and far field. In particular it very clearly brings out the similarity between the cohesive 
zone processes for opening or closing cracks while the contrasted effective adhesion 
energies are fully ascribed to the asymmetry found in the coupling. 

We believe these results pave the way to more extensive experimental comparisons 
between bulk viscoelastic response and effective adhesion of elastomers. 



Appendix A. Three-element model — Creep function 

The three-element model is an exponential creep function 

m = ^ (1 - (1 - A;) exp (-t/T)) (A.l) 

For the three-element model, 4>o(t), (fi^opif) and 0i,d(^) can be calculated 
analytically from Eqs. [H] [H] and [H] giving 

0o(t) = (l - ^ {1 - exp (-t/T)}j (A.2) 
0i,d(t) = (l - 2^/^ {1 - (1 + t/T) exp (-t/T)}^ (A.3) 
0i,op(t) = ^£7(1 + 2^y^ {1 - t/T - exp (-t/T)} j (A.4) 

Appendix B. Effective adhesion and effective dwell time 

It has been proposed [H], [13 [M] that the effective adhesion can be expressed as 

(receding contact) (B-1) 



Geff 



W Eoo^(j){t* 
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and 



— — = (growing contact) {^■^) 

where and 0(^6*) are effective compliances for the crack tip creep process during 
contact edge propagation, for opening and closing cracks. 

A linear theory of the coupling between local crack tip stress field and far field 
should involve the crack tip process in a similar way, irrespective of the contact edge 
direction (receding or growing). The fact that the effective crack tip compliance appears 
either at the numerator (Eq. lB.ll ) or the denominator (Eg. IB.2I) when the crack direction 
is reversed is therefore puzzling. 

Our equations Eqs. [19] and [2T] display a somewhat simpler structure. They explicitly 
call in the coupling (poitr) which results from the matching of the local stress field and 
the far field. Compared to Eq. IB.2t this form restores the effective crack tip compliance 
at the denominator of Eq. [21]so that the symmetry is preserved. In this respect, note that 
the Schapery Equations IB.2I and IB. II were first derived under the assumption of a power 
law creep function and presumably better apply to that case, where these equations 
clearly become equivalent to Eqs. [19] and [21] 
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